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Light top squarks (stops) and light higgsinos are the key features of natural SUSY, where the 
higgsinos xf and Xi 2 are nearly degenerate and act as the missing transverse energy {$rp) at the 
LHC. Besides the strong production, the stop can be produced via the electroweak interaction. The 
determination of the electroweak properties of the stop is an essential task for the LHC and future 
colliders. So in this paper, we investigate the single stop (ti) production pyp —>■ ti + $j< in the natural 
SUSY at the LHC, which gives the monotop signature t + $rp from fi —>■ txi ,2 or the monobottom 
signature b + from ii —>■ bx^. We perform Monte Carlo simulations for these signatures and 
obtain the results: (1) The signal b + $rj. has a better sensitivity than t + for probing natural 
SUSY; (2) The parameter region with a higgsino mass 100 GeV< < 225 GeV and stop mass < 

620 GeV, can be probed through the single stop production with S/'/B > 3 and 4% ^ S/B < 19% 
at 14 TeV HL-LHC with an integrated luminosity of 3000 fb“^. 


INTRODUCTION 

The search for supersymmetry (SUSY) is a long¬ 
standing important task in particle physics. One prime 
motivation for weak-scale SUSY is that it protects the 
Higgs vacuum expectation value without unnatural fine- 
tuning of the theory parameters. In the minimal super- 
symmetric standard model (MSSM), only a small subset 
of the supersymmetric partners strongly relates with the 
naturalness of the Higgs potential [1]. This can be seen 
from the minimization of the Higgs potential [2]: 

M| _ -P Ed) - {mjj^ -P S^) tan^ (3 2 

2 tan^ /3 — 1 ^ 
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where fx is the higgsino mass parameter, and and 
denote the weak scale soft SUSY breaking masses of 
the Higgs fields. A moderate or large tan/3 = Vu/vd is as¬ 
sumed in the last approximate equality. E^i and Tju arise 
from the radiative corrections to the Higgs potential, and 
the one-loop dominant contribution to E^ is given by [3] 

( 'm'x \ 

In order to obtain the observed value of Mz without large 
cancelations in Eq. (1), each term on the right hand side 
should be comparable in magnitude. Thus, the higgsino 
mass fi must be of the order of ~ 100 — 200 GeV and the 
requirement of E„ ~ Af|/2 produces an upper bound 
on the stop mass, which is about 500 GeV [4, 5] (a 125 
GeV Higgs mass can be achieved by a large stop tri- 
linear coupling without very heavy stops in the MSSM 


or achieved by extending the MSSM with additional D- 
terms or F-terms [6-8]). In addition, since the gluino 
contributes to at two-loop level, it is also upper 

bounded by the naturalness [9] (however, the direct LHG 
searches have pushed the gluino up to TeV scale [10, 11] 
while the recent ATLAS Z-peaked excess may indicate a 
gluino around 800 GeV [12]). 

So to test SUSY naturalness, the crucial task is to 
search for light stops or higgsinos. The search strat¬ 
egy for the pair productions of these nearly degenerate 
higgsinos has been recently studied [13-26]. During the 
LHC run-1, the ATLAS and CMS collaborations have 
performed the extensive searches for the stops through 
the gluino-mediated stop production [27, 28] or the direct 
stop pair production [29, 30]. Meanwhile, many theoret¬ 
ical studies that aim for improving the LHC sensitivity 
to a light stop have been proposed [31-46]. The current 
LHC constraints indicate a stop mass bound of hundreds 
of GeV [47-57], however, those results are affected by 
the stop polarization states and branching ratios. The 
constraints on the right-handed stop from the LHC run- 
1 direct searches [29, 30] are usually weakened by the 
branching ratio suppression, which can still be as light 
as 230 GeV in some compressed region [58]. If the stop 
mass is heavier than about 450 GeV, it is allowed in most 
parameter space. So in our work, we require the stop 
mass be heavier than 450 GeV [58]. 

Usually, the stop pair production provides the most 
sensitive way to search for the stop at the LHC. However, 
the stop can also participate in the electroweak inter¬ 
action processes. The determination of the electroweak 
properties of the stop is an essential task for the LHC 
and future colliders. In this work, we study the single 
stop production pp ^ ii in the natural SUSY at 
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the LHC. The observation of the single stop production 
will test the electroweak properties of the stop and the 
naturalness of the supersymmetry. In the following, we 
will perform the Monte Carlo simulations for the single 
stop production and examine its sensitivity at the LHC. 



Figure 1: The cross sections of t\t\ and tiyj” productions at 
the 14 TeV LHC for tan /3 = 50 and degenerate higgsinos with 
mass n = 100 GeV. The contribution of conjugate process 
i\xt included. 


CALCULATIONS AND SIMULATIONS 

At the LHC, the single stop production is induced by 
the electroweak interaction and proceeds through the fol¬ 
lowing process (see Fig. 1 for the corresponding Feynman 
diagrams): 

PP^iiXi- ( 3 ) 

Since in natural SUSY the light higgsinos are nearly de¬ 
generate, the decay products of xj" — W*Xi will carry 
small energies and, hence, are too soft to be observed in 
the detector. Thus, the associated production of tixT 
can be identified as ti + (!ti which provides a distinctive 
signature at the LHC. 

In Fig. 1, we show the next-to-leading order (NLO) 
cross sections of the stop pair and the single stop pro¬ 
ductions for /i = 100 GeV at 14 TeV LHC by using the 
packages Prospino2 [59] and MadGolem [60], respectively. 
The renormalization and factorization scales are taken as 
the half average of the final states masses. In the calcu¬ 
lations of tiXi , we use the LO and NLO parton densities 
given by CTEQ6L1 and CTEQ 6 M with five active flavors 
[61]. The contribution of the conjugate process i\xt is 
included. Except for the higgsino mass parameter /i and 
right-handed stop soft mass mu ^, we assume other soft 
supersymmetric masses at 1 TeV, and use the packages 
SOFTSUSY-3.3.9 [62] and MSSMCalc [63] to calculate 
masses, couplings and branching ratios of the sparticles. 


Since the cross section of single stop production is not 
sensitive to tan ,5, we take tan ,5 = 50 for simplicity. We 
find that the single stop cross section can still reach about 
200 fb when ~ 450 GeV. The NLO AT-factor of the 
process pp —>■ Uxr ranges from 1.25 to 1.33. When the 
stop becomes heavy, the single stop production cross sec¬ 
tion will decrease, but slower than the pair production, 
due to the kinematics. 

Next, we investigate the LHC observability of the sin¬ 
gle stop signatures with the sequent decays ti txi 2 
and ii —>■ 6 x 1 " • 

PP iiXi tXp2Xi bjj + (4) 

PP^hXx ^^xtX\ +t'T- (5) 

For the decay ii tXi, 2 i the SM backgrounds to the 
signal bjj + S't £^re from the semi- and full-hadronic ti 
events [64-66] ^, where the undetected lepton and the lim¬ 
ited jet energy resolution will lead to the relatively large 
missing transverse energy. The processes W + jets and 
Z -|-jets can fake the signal when one of those light-flavor 
jets are mis-tagged as a 6 -jet. Also, the single top can 
mimic our signal when the lepton from the W boson de¬ 
cay is missed at the detector. While ti +V backgrounds 
are not considered in our simulations due to their small 
missing energy or cross sections compared to the above 
backgrounds. 



p^{top)/GeV 


Figure 2: The parton-level pr distribution of the top quark 
in the channel ti tXi ,2 for p = 100 GeV at 14 TeV LHC. 


^ Hadronic monotop has received special attention since its signa¬ 
ture offers the possibility to use top reconstruction as a tool to 
reject the backgrounds. In contrast, the leptonic monotops is 
believed to be more challenging since the branching fraction of 
the leptonic top quark decay is smaller and since there are two 
different sources of missing transverse energy, namely a neutrino 
coming from the top quark decay and the new invisible state. In 
Ref. [65], the authors comparatively studies these two channels 
and found that the sensitivity of both channels are very similar. 
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In Fig. 2, we present the parton-level px distribution 
of the top quark in the channel ti —>■ 2 = 100 

GeV at 14 TeV LHC. It can be seen that, with the in¬ 
crease of stop mass, the top quark produced from stop 
decay is boosted and has larger px- So, in the analysis of 
ii 2 channel, we respectively adopt HEPTopTagger 

[32] and normal hadronic top reconstruction methods for 
each sample to identify the top quark in the final states 
and present our results with the best one. The detailed 
analysis strategies are the followings: 

• Events with any isolated leptons are rejected; 

• Method-1: We use Cambridge-Aachen (CA) algo¬ 
rithms [67] in Fastjet [ 68 ] to cluster the jets with 
R = 1.5 to obtain the top-jet candidates. Each 
candidate must have the top quark substructure 
required by the HEPTopTagger. The 6 -tagging is 
also imposed in the top-jet reconstruction. Other 
energy deposits outside the top-jet are further re¬ 
constructed as the normal jets by using anti-fct al¬ 
gorithm with R = 0.4 [69]. The top window used 
in our analysis is 150 < mt < 200 GeV. While the 
W window is taken as the default value in HepTop- 
Tagger; 

• Method-2: In normal hadronic top quark recon¬ 
struction, a pair of jets is selected with the invariant 
mass rrijj > 60 GeV and the smallest AR. A third 
jet closest to this di-jet system is used to consti¬ 
tute the top quark candidate. Among these three 
jets, at least one 6 -jet and A(j){$xiPT{bi)) > 1 is 
required. The anti-fc* algorithm is used for jet clus¬ 
tering with i? = 0.4; 

• We keep the events with the exact one recon¬ 
structed top quark and require 150 GeV < < 

200 GeV; 

• The extra leading jet ji outside the reconstructed 
top quark object is vetoed if px{ji) > 30 GeV and 
l^0-i)l <2.5; 

• We define the signal regions according to 
{tx,PTijtop)) cuts: (200, 100), (250, 150), (300, 
200), (350, 250) for Method-1, and {px{b),$x) 
cuts: (200, 50), (250, 50), (300, 100), (350, 100) 
GeV for Method-2. 

For the decay H —>■ , the SM backgrounds to the 

signal 6 -F $x are dominated by the processes W -F jets 
and Z -F jets when the light-flavor jets are mis-identified 
as 6 -jets [70]. The tt events become the sub-leading back¬ 
grounds due to their large multiplicity. The signal events 
are selected to satisfy the following criteria: 

• Events with any isolated leptons are rejected; 


• We require exact one hard 6 -jet in the hnal states, 
but allow an additional softer jet with px{ji) < 30 
GeV and A(j){$XTPT{ji)) > 2. Since the hardness 
of 6 -jet from stop decay depends on the mass split¬ 
ting between ti and x]", we define three signal re¬ 
gions for each sample according to {$x^PT{b)) cuts: 
(100, 70), (150, 100) and (250, 200) GeV. 

Finally, we use the most sensitive signal region (with the 
highest S/y/B) for each decay mode and show our re¬ 
sults in Fig. 3 and Fig. 4, respectively. In our study, we 
omitted the QCD multijet backgrounds, whose correct 
treatment needs the experimental data-driven methods 
and hence depends on the realistic detector environments 
of the 14 TeV LHC. As discussed in [65, 70, 71], the re¬ 
quirements of high pt ( 61 ) and large $x with a separation 
Acl){$xiPT{ji)) can usually be expected to greatly re¬ 
duce the fake contamination from the QCD backgrounds 
and allow for a good signal selection efficiency. 

The parton level signal and background events are gen¬ 
erated with MadGraph5 [72], where W/Z-Fjets is matched 
up to 3 jets by using MLM matching scheme [73] and set¬ 
ting xqcut = 30 GeV. For the value of qcut in matching, 
we take it to max{xqcut-\-5, xqcut *1.2) [74] in our simu¬ 
lation. We perform parton shower and fast detector sim¬ 
ulations with PYTHIA [75] and Delphes [76]. We assume 
the 6 -jet tagging efficiency as 70% [77] and a misiden- 
tification efficiency of c-jets and light jets as 10 % and 
0.1%, respectively. The cross section of ti is normalized 
to the approximately next-to-next-to-leading order value 
au = 920 pb [78]. 

RESULTS AND DISCUSSIONS 


Table I: The cross sections of V-Fjets, U and ti(— tXi, 2 )Xi for 
a benchmark point = (611,100) GeV and tan/1 = 10 

in Method-1 and Method-2 at 14 TeV LHC with C = 3000 
fb“^. The cross sections are in unit of fb. 


cuts 

W +jets 

Z + jets 

ti 

tw 

S 

SjB 

s/Vs 

Method-1 

< 10"'' 

0.29 

2.20 

0.80 

0.13 

4.0% 

3.9 

Method-2 

< 10"'' 

0.59 

0.55 

0.24 

0.044 

3.2% 

2.1 


In Table I, we compare the results of pp —> ti(—>■ 
tXi, 2 )xr for a benchmark point (m^^,/*) = (611,100) 
GeV in Method-1 and Method-2 at 14 TeV LHC. From 
this table we can see that the Z -F jets background in 
Method-1 is smaller than in Method-2, while the ti back¬ 
ground in Method-1 is larger than in Method-2. However, 
the signal events can be more kept in Method-1 than in 
Method-2. So the overall effects make the Method-1 have 
a better sensitivity in reconstructing the top quark in the 
region with large mass splitting between ti and Xi- At 
14 TeV LHC with C = 3000 fb“^, the statistical signif- 
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icance Sj'/B for our benchmark point can reach 3.9cr 
(2.1cr) with S/B = 4.0% (3.2%) in Method-l(2). 

Table II; The cross sections of fo+jets, ti and ti(—>■ bxt)Xi for 
a benchmark point — (496, 200) GeV and tan/3 = 10 

at 14 TeV LHC with £ = 3000 fb~^. The cross sections are 
in unit of fb. 


W +jets 

Z -b jets 

tt 

S 

S/B 

S/y/B 

< 10“'' 

2.77 

1.10 

0.20 

5.1% 

5.5 


In Table II, we show the cross sections oi V + jets, 
ti and ti(—f bxt)Xi for a benchmark point (m^^,//) = 
(496, 200) GeV at the 14 TeV LHC. Different from ii —>■ 
txi 2 channel, Z + jets background is dominant over ti 
since only one hard 6-jet is required in the final state. 
From Table II we can see that S/y/B and S/B can reach 
5.5 and 5.1% for £ = 3000 fb“^, respectively. 



500 550 600 650 700 750 800 


riTj (GeV) 

Figure 3: The dependence of the significance of the channel 
ii —>■ tXi ,2 on the higgsino mass /i and stop mass mi_^ at the 
14 TeV LHC with £ = 3000 fb"b 

In Fig. 3, we display the dependence of statistical sig¬ 
nificance S/yfB of the channel £ —?► txi 2 on the higgsino 
mass n and stop mass at 14 TeV LHC with £ = 3000 
fb“^. We can see that values of S/yfB decrease with 
the increase of ^ because of the cut efficiency reduction. 
When the stop becomes heavy, the cross section of Cxj" 
is suppressed. However, as a result of the application of 
HEPTopTagger method, more signal events can be kept, 
in particular in the mass range 450 GeV < rrif^ < 650 
GeV. Therefore, when ^ < 150 GeV, the stop mass < 
610 GeV can be probed at > 3a statistical significance 
with S/B < 8%. 

In Fig. 4, the statistical significance S/y/B of the chan¬ 
nel £ —>■ bx// is presented on the plane of higgsino mass 
/r versus stop mass at 14 TeV LHC with £ = 3000 
fb“^. It can be seen that the sensitive stop region lies in 
450 GeV < < 620 GeV, where a hard 6-jet {pr > 200 

GeV) and the sizable > 250 GeV) can be used to 


LHC-14 TeV, Lum.=3000 fb ’ S/Vb 



trij (GeV) 

Figure 4: Same as Fig. 3, but for the decay channel H —>■ bx/ ■ 

effectively suppress the backgrounds. But when the stop 
mass increases, S/yfB will rapidly decrease. We see that 
the higgsino mass 100 GeV< p < 225 GeV and the stop 
mass < 620 GeV can be covered at > 3a statistical 
significance with S/B varying from 4% to 19%. 

CONCLUSIONS 

even its exclusion limit for the stop may be not as 
good as the pair production. However, the stop can par¬ 
ticipates in the strong interaction processes but also in 
the electroweak interaction processes, the determination 
of the electroweak properties of the stop is an essential 
task for the LHC and future colliders. In this work we 
propose to probe natural SUSY by using the electroweak 
single-stop productionpp —>■ ii+^rp at the LHC (here the 
missing energy is from the nearly degenerate higgsinos). 
By analyzing the decay channels of the stop ii —>■ 6%° 2 
and ii —>■ bx//, we obtain the obervations: (1) The de¬ 
cay ii -X bx/ has a better sensitivity than H ty? 2 
; (2) The parameter region with a higgsino mass 100 
GeV< pL < 225 GeV and the stop mass < 620 GeV 
can be covered with S/y/B > 3 and 4% ^ S/B < 19% at 
14 TeV HL-LHC with an integrated luminosity of 3000 
fb“^. So the searches for the single stop production will 
directly test the naturalness of the supersymmetry and 
the electroweak properties of the stop. 
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